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ABSTRACT 

A study is made of the ionization of the Martian atmosphere below 100 km 

i 

in the context of recent data on the atmosphere and ionosphere of the planet. 
Suggested sources of ionization include: solar ultraviolet and x-radiation, 
cosmic rays, meteors, radioactive material within the atmosphere and radio- 
activity of the surface, the solar wind, and protons originating from solar 
flares. Loss processes, including charge rearrangement, dissociative re- 
combination of ions and electrons, and negative ion formation are detailed. 
Possible importance of the negative ion CO^ is discussed. Its distribution 
will depend on the0 2 content of the atmosphere as well as the distribution of 
O and CO. CO - and CO ~ may also be present in the ionosphere. It is shown 
that combination of the solar photoionization rate with that due to cosmic radia- 
tion will lead to an ion density minimum. The magnitude of the ion density and 
location of the minimum will depend on the atmospheric model chosen, solar 
activity, and solar zenith angle as well as such factors as the presence of par- 
ticulate matter and meteor ionization. 


3 


THE LOWER IONOSPHERE OF MARS 


INTRODUCTION 

Recent experimental results have yielded much definite information on the 
properties of the Martian atmosphere. In particular spectroscopic studies have 
given us a more complete picture of the chemical composition, temperature, 
and pressure distribution with altitude, Gray (1966), Belton and Hunten (1966). 
Similar information has been derived from the Mariner IV occultation experi- 
ment conducted during the fly-by of Mars in July 1965, Kliore et al (1965). 
Mariner IV also yielded important data concerning the ionosphere of Mars, 
Fjeldbo et al (1966a,b), Kliore et al (1965). The electron concentration distri- 
bution as derived from the Mariner IV 2 frequency occultation experiment is 
shown in Figure 1, Fjeldbo et al (1966a). An electron density scale height of 
24 ± 3 km was found together with a maximum electron density of 9.5 ± 1.0 x 
10 4 el/cm 3 at an altitude of 120 + 5 km. The plasma height remained constant 
to an altitude of 250 km. There is a minor ionospheric layer at 100 km with a 
maximum electron density of 2.5 x 10 4 cm’ 3 . An upper limit of 5 x 10 3 cm -3 
was found for the nighttime ionosphere, Kliore et al (1965), Fjeldbo et al (1966a). 

Several atmospheric models have been proposed to explain the occultation 
data. These can be categorized according to their description of the ionosphere 
as analagous to the Earth's ionospheric E, F-l, and F-2 regions as regards to 
the formation. The designations are interpreted to mean that in the E -region 
case the ionospheric peak is formed by a balance the ion pair production rate 
due to long wavelength solar x-rrays with recombination of CO* ions, Chamberlain 
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and McElroy (1966) , Me Elroy (1967) . The primary ion is considered to be 0 + 
at the observed electron density maximum of 120 ± 5 km for the F-l region, 
Donahue (1966). Diffusion of ions and electrons to form the peak is the basis 
of the F-2 region, Johnson (1965). 

The present paper is a study of the lower ionosphere of Mars, i.e., that 
portion which lies between the surface of the planet and an altitude of approxi- 
mately 100 km. The description of such an ionosphere will depend on the 
atmospheric model assumed. For comparison, two atmospheric models are 
employed, namely, the Johnson "F-2" and the Chamberlain and McElroy "E". 

It will be shown that an F-l type atmosphere gives rise to a lower ionosphere 
between these two extremes. 

SOURCES OF IONIZATION 

The following sources of ionization must be considered for the atmosphere 
of Mars: 

1. Solar ultraviolet and x-radiation 

2. Cosmic rays 

3. Meteors 

4. Natural radioactivity of the planetary surface and radioactive material 
within the atmosphere. 

5 . The solar wind 

6. Protons originating from solar flares 

Since the planetary magnetic field is at most 3 x 10~ 4 times that of the 
Earth, Smith et al (1965), no radiation belts and auroral particle precipitation 
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zones are to be expected. Indeed Van Allen et al (1965) and O' Gallagher and 
Simpson (1965) found no evidence of trapped radiation. 

In considering the effect of solar radiation it is essential to know the in- 
tensity as a function of wavelength, as well as the cross sections for absorption 
and photoionization of the gases constituting the Martian atmosphere. The in- 
tensity of the solar radiation at 1AU has been summarized by Hintereger, et al 
(1965). The intensity at distances other than 1AU can easily be calculated 
knowing the properties of the orbit of Mars. The distance of the planet from 
the sun varies between 1.36 and 1.67 AU over a two year period. The mean 
distance is 1.52 AU. This leads to a dilution factor of between 0.524 and 0.36, 
Schultz and Holland (1962). At the time of Mariner IV the planet was in a posi- 
tion such that the dilution factor was 0.39. 

As stated previously the possible constitutents of the Martian atmosphere 
are C0 2 , N 2 , Ar,0,C0 and possibly0 2 and oxides of nitrogen such as NO. 

The wavelength equivalent of the ionization potential of these constitutents is 
given in Table I. 


Table I 

Ionization Thresholds of Martian Atmospheric Constituents 


Ar 

787A 

n 2 

7 96 A 

CO 

885A 

co 2 

899A 

o 

91 0A 

o 2 

1026A 

NO 

1340A 
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can be overcome if an additional ionization source exists. One source not 
previously considered is that of meteor ionization. In the Earth's atmosphere 
most meteors deposit their energy between 90 and 120 km with some coming 
as low as 70 km. For an incoming meteor of mass m, density p mt and velocity V, 
the number of ion-pairs produced per unit path length is according to McKinley 
(1961) 


q 


m 

r B — 
q A 


2/3 


Pa V 4 


m 


where r q is the ionization efficiency factor B is a constant which is a function 
of the shape and ablation properties of the meteor and p m is the density of the 
atmosphere. Assuming that the atmospheric density distribution is described 
by a constant scale height H the maximum ionization will occur where 


p a H - C m 1/3 P m 2/3 V" 2 cos X 


Here the path length (ds) in terms of the verical height, dz, is 

ds = dz/cos X 


and C is another constant which is a function of B. Table II gives the altitude 
of maximum ionization for the different Martian atmosphere models in com- 
parison with deposition altitudes for equivalent meteors impinging on the Earth's 
atmosphere. 
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Some of the most intense line and continuum emissions of the solar ultra- 
violet spectrum occur at wavelengths longer than 900A. This radiation includes 
the important Lyman series Lyman alpha (1215.7), Lyman beta (1025.7), etc. 
and the Lyman continuum (910A) which will be almost ineffective in producing 
ionization unless NO,0 2 , and O are present in the atmosphere. For an atmosphere 
composed principally of C0 2 , Hel (584. 3A) and Hell (303. 8A) as well as con- 
tinuum radiation below 200 A will play an important role. 

Figure 2 illustrates the ion-pair production functions due to solar radiation 
for the McElroy E and F-l models and the Johnson F-2 model for a solar zenith 
angle of 67°. The latter was calculated by the author; the E and F-l models are 
reproduced from McElroy (1967). The peak production rates correspond in all 
cases to unit optical depth for wavelengths in the vicinity of 304A. The second- 
ary peaks at lower altitudes correspond in the E and F-l region cases long 
wavelength x-rays. For the F-l model the peak of production should corre- 
spond to the peak of the electron density distribution provided the loss processes 
are not altitude dependent. McElroy (1967) has demonstrated that in the case 
of the F-l layer, the maximum of ionization occurs below the production peak. 

The same thing occurs in the case of the E region, the 120 km peak being 
caused by solar radiation at x-ray wavelengths, where the flux is extremely 
variable. The F-2 model has the peak of production near 100 km. Such a 
model requires a loss rate in excess of 10" 6 cm 3 sec -1 to explain the observed 
low altitude variation. 

The problem of the need for an additional ionization source at low altitudes 
which is most acute for the E region model, but exists for the F-l model also, 
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dust should be at least as probable on Mars as on the Earth. The actual pres- 
ence of meteoric layers (already established for the Earth's atmosphere) will 
depend on the loss process acting in the Martian atmosphere. 

In addition to direct ionization by meteor ablation, photoionization of accumu- 
lated meteoric atoms can occur. These atoms are primarily Fe, Mg, Si, Ca, 
and Na. Their presence in ionized form has been detected in the Earth's atmos- 
phere at the turbopause by rocket-borne ion mass spectrometers both during 
ordinary and meteor shower conditions. Narcisi et al (1967 a,b). However, 
estimates of the ionization rate coefficients were given many years previously 
by Nicolet (1955). These rates I(Na) = 10 -s sec -1 , 1(Ca) = 10 -6 sec -1 , 1(Mg) = 

10 -7 sec -1 can be used to estimate the number of atoms necessary in order to 
yield an ion density of 10 3 cm -3 , assuming a radiative recombination coefficient 
of 10 -12 cm 3 sec -1 . The density of atoms is less than 10 cm" 3 . Fehsenfeld et al 
(1968) have given rate coefficients for the following processes 


Mg + + 0 2 + Ar 

—MgO, 

+ Ar 

Ca + + 0 2 + Ar 

,CaQ 2 

+ Ar 

Fe + + O a + Ar 

_*Fe0 2 

+ Ar 

Na + + O 2 + Ar 

,Na0 2 

+ Ar 


Similar reactions may occur in the case where CO a is substituted for 0 2 . Such 
reactions or an atmosphere containing sufficient 0 2 could lead to molecular 
ions. These ions will dissociatively recombine with electrons with a rate coeffi- 
cient in the 10 -6 to 10 -7 cm 3 sec -1 range. Under such circumstances metalic 
atom densities in the 10 6 to 10 7 cm * 3 range are necessary in order to produce 
10 3 electrons cm -3 . 
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Table H 


Deposition Height of Meteors 


Earth Altitude 

McElroy E 

McElroy F-l 

Johnson F-2 

70 

74 

61 

44 

80 

91 

74 

55 

90 

111 

89 

65 

100 

129 

102 

73 

110 

145 

113 

80 


Various estimates of the contribution of meteors to the ion-pair production 
function have been made for the earth and range in value between 10“ 5 and 10" 1 
electrons cm -3 sec" 1 e.g. Lovell (1957). Fiocco (1967) has reexamined the 
problem considering recent information on the interaction of hypersonicly 
moving bodies with the Earth's atmosphere. The ionization rate depends on 
the area/mass ratio of the particle as well as its velocity, the efficiency of 
ionization and flux of incoming particles. For velocities between 39.2 and 44.7 
km/sec and area/mass ratios between 10 and 10 3 m 2 /kgm, the production 
function varies between 10 15 and 10 18 ion pairs/kgm m with values between 
10 17 and 10 18 /kgm m being most typical. 

An average value for the incoming flux at the Earth is 10~ 13 kgm/m 2 sec, 

i * 

which could lead to an ion-pair production function of 10“ s /cm 3 . Dust measure- 
ments by Alexander et al (1965) conducted on Mariner IV indicated a factor of 
five increase in the flux at 1.4 AU as compared with that at 1 AU, although there 
was a slight decrease in the immediate vicinity of Mars. Thus ionization by 
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An ionization chamber was included as part of the instrumentation of 
Mariner IV. The instrument ceased operation on 16 February 1965. One can 
compare this data with that obtained by the Mariner II probe to Venus in 1962. 
The data, taking into account the masking effect of the planet, gives 330 ion 
pairs/cm 3 sec atm for 1962, Neher and Anderson (1964), and 500 ion pairs/cm 3 
sec atm for 1965, Neher as reported by Nicks (1967). Assuming the ionization 
efficiency of C0 2 to be the same as air, the ionization rate for the time of 
Mariner IV can be written as 

q CR (z) = 1.85 x 10" 17 n(M) ion pairs/cm 3 sec 

where n(M) is the local number density at any altitude, Z. Figure 2 illustrates 
the variation of q c r as a function of altitude for the different atmospheric 
models. 

The undisturbed solar wind is composed of protons of velocity averaging 
between 300 and 400 km/sec (500 to 900 ev) with densities up to 20/cc (Parker, 
1964). Because of the absence of a magnetic dipole of strength greater than 
10“ 3 of the Earth; the solar wind will penetrate the Martian atmosphere. Indeed 
Sagan and Veverka (1967) explain the subsidiary ionospheric peak as being caused 
by the solar wind impinging on an F-2 type atmosphere. Recent information 
from the Mariner V probe to Venus give strong evidence that the solar wind 
does not penetrate the atmosphere. These results indicate an ionospheric 
conductivity-solar wind magnetic field interaction which causes the solar wind 
to slip around the planet. The situation should be similar for Mars, Mariner 
Stanford Group, (1967), Lazarus el; al (1967). 
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Table in illustrates the altitude to which more energetic protons will pene- 
trate for different Martian atmospheric models. Solar flare protons will dissi- 
pate most of their energy near the surface of the planet. 

Table HI 


Depth of Penetration of Protons into the Martian Atmosphere 


Model 

10 Kev 

100 Kev 

1 Mev 

10 Mev 

50 Mev 

Johnson F-2 

80 km 

72 km 

57 km 

35 km 

11 km 

McElroy F-l 

103 km 

88 km 

70 km 

35 km 

10 km 

McElroy E 

140 km 

125 km 

95 km 

50 km 

15 km 


As a result of the natural abundance of the elements , it would be expected 
that the ion-pair production due to the decay of elements within the Martian 
surface should be comparable to that of the Earth. It has been noted previously 
that the main production of cosmic ray secondaries is below the Martian surface. 
This may increase the radioactivity of the surface. However, for our purposes 
we shall assume that the ion-pair production function is the same as that of the 
Earth. Bricard (1965) has given an excellent review of the role of radioactivity 
in atmospheric electricity and estimates that the contribution of the surface i s 
3.5 + 2.5 - 1.5 ion-pairs/cm 3 sec/atmosphere. 

According to Bricard (1965) (Table I after Israel (1961)), 0.1% of the soil 
radiation remains at 1 km above the surface of the Earth. This corresponds to 
a range in excess of 100 gm/cm 2 , surface radiation will suffer little attenuation 

in its passage through the atmosphere of Mars. At the same time it represents 

\ 

only about 1% that which is produced by cosmic radiation. At this point it is not 
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possible to speculate on the radioactivity content of the Martian atmosphere. 
Radioactivity aerosols may well exist. 

POSITIVE ION REACTIONS 

The principal positive ion reactions \yhich will occur in the Martian iono- 
sphere are listed in Table IV together with their rate coefficients and sources. 
Also listed are the important positive ion- electron reactions which should be 
considered. For an atmosphere where 0 + is the predominant ion produced by 
photoionization C0 2 will result from the reaction of 0 + with C0 2 . Rapid charge 
exchange will occur between N 2 and C0 2 . Since the dissociative recombination 
coefficient of C0 2 as measured by Weller and Biondi (1967) is 3.8 ± 0.5 x 10~ 7 
cm 3 sec" 1 at 300°K and of the same order at 210°K, 0 2 concentrations in excess 
of 4 x 10 8 cm" 3 must be present for electron densities of 10 s cm” 3 or greater 
in order for 0 2 to become the dominant ion when reactions between C0 2 and 0 2 
occur. 

Table IV 

Important Positive Ion Reactions in the Martian Ionosphere 

Reaction Rate 

cm 3 sec" 1 

Ar + + C0 2 — ►Ar + C0 2 d) 7.4 x lo" 10 C + + C0 2 -^C0 + CO + (8) 1.8 x 10' 9 

Ar + + CO— Ar + CO + (l) 8 .OXIO " 11 C + +0 2 — -O + CO + (8) l.l x io“ 9 

Ar + + 0 2 — -Ar + O* (i) 1.1 xio -10 

Ar + + N 2 — ►Ar + N 2 (i) 6 x 10 " 12 

C0 + +C0 2 — C0+C0 2 (8) 1.1X10" 9 

C0 + +0 2 ►CO+Oj (8) 2X10" 10 



Table IV (Continued) 

Important Positive Ion Reactions in the Martian Ionosphere 
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Reaction 


Rate 

Nj + C0 2 — ►Nj + C0 2 


8 x 10" 10 

n; + o 2 - 

— n 2 + o + 2 

(2) 

1 x lO -10 

N + 2 + O — 

— ►N + NO + 

(3) 

2.5 x 10' 10 

0*2 + N 2 - 

— ►NO + + NO 

(2) 

>10' 15 

C0 2 + 0 2 — ►COj + o 2 


1.0 x 10' 10 

C0 2 + NO 

-►C0 2 + NO + 


— 

o + + co 2 — ►CO + o 2 

(8) 

1.2 x 10~ 9 

o + + o 2 — 

— ►Oj + O 

(2) 

2 x 10" 11 

o + + n 2 — 

— ►N + NO + 

(2) 

4.6 x 10“ 12 

C0 2 + e — ►CO + O 

(4) 

4 x 10“ 7 

0 2 + e — 

—►O + 0 


5.2 x 10' 6 /T e ° 

N 2 + e — 

-►N + N 

(5) 

2 x 10“ 6 /T e °- 

NO + + e - 

— ►N + O 

(6) 

1.5 x 10" 4 /T 

CO + + e — 

-►c + O 

(7) 

6.8 x 10“ 7 

Ar + + e — 

— ►Ar + hv 


1 x 10' 12 

0 + + e 

-►O + hv 


1 x 10“ 12 

C + + e 

— ►C + hv 


1 x 10~ 12 

Mg + + e 

“►Mg + hv 


1 x io' 12 


1) Warneek (1967a); 2) Warneck (1967b); 3) Ferguson et al (1965); 

4) Weller and Biondi (1967); 5) Mehr and Biondi (1968); 6) Weller and 
Biondi (1968); 7) Mentzoni and Donohoe (1968); 8) Norton et al (1966). 
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The measured rate of the C0 2 dissociative recombination coefficient, 

4 x 10" 7 cm 3 sec” 1 does not support the assumption of 10" 5 cm 3 sec” 1 made by 
McElroy (1967) in order to obtain agreement between the E region model and 
the experimental data. The assumption of a positive exponent temperature 
dependence of a 0 (10” 8 cm 3 sec” 1 at 330°K and 10” 7 at 300°K) is also in dis- 
agreement with laboratory data. 

In contrast the F-2 model of Johnson requires that the a Q be of the order 
of 10' 5 cm” 3 sec” 1 where the E region would have it be 10“ 8 cm” 3 sec” 1 . In 
the present work it will be assumed that a Q ^ follows the temperature dependence 
given in Table IV and that a c o 2 has a T” 0 - 5 orT” 1,2 dependence. This choice 
is consistent with laboratory measurements. 

NEGATIVE ION REACTIONS 

Laboratory experiments have not found direct attachment of electrons to 
C0 2 at thermal energies, Pack et al (1962) and Schlumbohm (1962). Attachment 
has been found for electron energies in excess of 2 ev, Schulz (1962). Ion mass 
spectrometer analysis has been used by Moruzzi and Phelps (1966) to establish 
the attachment process at these ionosphericly large electron energies. The 
reactions 

e + C0 2 *CO + O 

and 

0”+2C0 2 ►C0 3 ” +C0 2 

were found to occur. The ion C0 3 is found to be dominant and would be ex- 
pected in a pure C0 2 atmosphere with low energy electrons. 
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Several other species are possibilities for trace amounts of constituents 
common to the atmosphere of the Earth. One such example is the formation of 
CO” which will occur by the three-body process involving thermal electrons 

A 

e + 2CO ►CO - + CO 

at a rate of 2.7 x 10~ 31 cm 3 sec -1 Mentzoni and Donohoe (1968). The reactions 
resulting in the loss of this ion are unknown at present. However, this ion is a 
possibility if CO is present in the lower atmosphere of Mars. 

An example of a negative ion distribution can be given for an atmosphere 
which contains 0 2 . The three-body attachment of thermal electrons to form 
0 2 - is well established, Chanin, Phelps, and Biondi (1962). This process 

e + 0 2 + M ►0 2 - + M 


is followed by 


0 2 +C0 2 +M ►C0 4 " + M 

Pack and Phelps (1966a) have determined that when M is C0 2 , the rate coefficient 
is 2 x Hr 30 for the attachment process. Mourzi ahd Phelps (1966) have deduced 
a rate of 9 x lO~ 30 cm 6 /sec for C0 4 formation where M is C0 2 . If M is 0 2 , 
this rate coefficient is 2 x 10 -29 cm 6 /sec. 

Other charge exchange reactions which can occur in 0 2 , C0 2 nurtures 


include 
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O +O — 

*-0 2 + e 

1.9 x l(f 

0 + o 2 — 

— + 0 

— 

O - + 2 0 2 - 

— o 3 - +o 2 

4 x 10 -31 

O3 + co 2 

-CO3- +o 2 

4 x lo -10 

co 3 - + 0 

— -o 2 + co 2 

8 x IQ -11 


For an atmosphere having O as a principal constitutent direct attachment 
takes place with a rate of 1 x 10 - 15 cm 3 /sec. 

Processes leading to the destruction of negative ions must include ion-ion 
recombination, associative detachment, collisional detachment and photo- 
detachment. The electron affinities of many of the molecules which result in 
negative ions are unknown at present. The affinity of O is reasonably estab- 
lished at 1.5 ev, Schulz (1962). A lower limit of 0.15 ± 0.05 ev was set by Burch 
et al (1958) in order to explain photodetachment cross section data of O J . 
Phelps and Pack (1966) obtained a value of 0.43 ± 0.02 ev. The affinity of CO^ 
is 1.22 ev Pack and Phelps (1966). The affinities of CO J and CO" are unknown. 
One would expect that visible light would not result in detachment of electrons 
from CO J and CO 4 . In addition to needing shorter wavelengths one must con- 
sider the reduced intensity which occurs at Mars. 

Associative detachment can be an important loss process particularly when 
O and CO are present. The reaction 

CO + O" -C0 2 + e 

has a rate of 8 x 10 - 10 cm 3 sec -1 , Moruzzi and Phelps (1966). In addition the 


reactions 
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and 


Oj + O— ►Oa+e 3.3 xio -10 


O + O— 0 2 + e 1.9 x io' 10 

must be considered. For rate Coefficients of this order negative ions will be 
destroyed for all but the smallest concentrations of O. 


Collisional detachment was not found to be appreciable in 0 2 — CO a mixtures 
where C0 2 is the major constitutent. However the reaction 

C0 4 " + C0 2 ►2C0 2 + 0 2 + e 

does occur at high pressures Mourzi and Phelps (1966). 


The rate coefficients for such ion - ion recombination processes as 


and 


C0 3 " + C0 2 + -2C0 2 + O 

►CO + co 2 + 0 2 


C0 4 ' + C0 2 + ^2002 + 0 2 


as well as similar negative ions reacting with C0 + and 0 2 are unknown at 
present. 


On the basis of the abundance of constituents of the Martian atmosphere, 
the following negative ion reactions may be selected as being the most im- 
portant 
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e + 0 2 + C0 2 — 

— o 2 + co 2 

a 

0 2 + hv 

— *“ 0 2 + e 

d 

0 2 + 2C0 2 

-►CO 4 + co 2 

k i 

C0 4 + hv 

— ^C 0 4 + e 

d i 

CO 4 + co 2 

— *-2C0 2 + Oj + e 

k 2 

CO 4 +C0 2 + 

— *-2C0 2 +0 2 

a. 

1 


Ion - ion recombination will become important only during disturbed events and 
can be neglected in the present analysis. The resulting equilibrium expression 
for the ratio of CO 4 ions to electrons, N e , is 

ak, [Q a ] [C0,]» 

C0 * {d + k! [C0 2 ]^} {d,*k, [CO,]} 

Molecular oxygen, while not having been detected spectroscopicly, is prob- 
ably a minor constituent in the atmosphere. Berkner and Marshall (1967) 
estimate that the upper limit is 1 % of the total atmosphere. The E region model 
assumes the concentration of 0 2 to be 10" 4 of the total number density. 

Figure 3 gives X . c0 - as a function of C0 2 number density for an assumed 0 2 
concentration of 10 ” 2 [C0 2 ]. A value of 2 x 10" 17 cm 3 sec -1 will be taken for 
k 2 and Q.2 sec -1 will b$ chosen for the photodetachment rate of O'. Photo- 
detachment from CO 4 will be neglected. 

The dashed curve corresponds to the nighttime case of no photodetachment. 
In such circumstances 



= 10 ' 15 co 2 


3 [° 2 ]/ k 2 
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The altitude for k = 1 can occur at about the same altitude as on Earth for the 
E region model. For an F-2 model \ = 1 almost 20 km below the corresponding 
E region altitude. In any event, for these assumptions k = 1 for [C0 2 ] = 

10 15 cm -3 . 

DISTRIBUTION OF IONS AND ELECTRONS 

The results of the previous sections can now be combined to estimate the 
ion density distribution with altitude employing the sources of ionization and 
loss processes discussed earlier. In Figure 4 is shown the lower Martian 
ionosphere as it would appear if the observed electron density distribution is 
that of an E region. The observed electron density distribution will decrease 
rapidly to a minimum of the order of 100 cm -3 at 85 km. Below this altitude 
the positive ion density increases due to cosmic ray ionization, the minimum 
represents the matching of the solar photoionization rate with that of the cosmic 
rays. Ionization increases to the planet's surface where the ion density is of 
the order of 3 x 10 3 cm -3 . Also shown is the negative ion, CO^ , density cor- 
responding to the curve of k given in Figure 3. Because of the reaction between 
C0 2 + and 0 2 the major ion of the cosmic ray layer will be O 2 even if 0 2 is only 
10" 4 of the total number density. 

Figure 4b shows the distribution which would be expected if the ionosphere 
is an F-2 layer. In this instance the minimum will occur in the vicinity of 50 km. 
The species of ion will be the same as in the E region case. The effect of 
meteors has not been included in these ion density estimates. In the case of 
the E region model the deposition heights are such that the meteors contribute 
mostly in the height range where the maximum of ionization is observed. In 
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the F-2 case the deposition will occur in the 75 to 85 km range — see 
Table II. 

PARTICULATE MATTER IN THE ATMOSPHERE OF MARS 

The existence of a particulate matter on Mars is well established. Opik 
(1966) has reviewed the optical observations of the planet and concludes — 
"Limb darkening at full phase in green, yellow, and red light indicates absorp- 
tion by atmospheric haze, aerosols, and dust." An average particle radius of 
10' 3 cm was deduced from the lifetime of the large (area 3 x 10 3 km) yellow 
cloud which was observed in 1956. Greenspan (1965) has postulated that con- 
centrations of 10 8 to 10 9 ice particles/ cm 2 column and of 0.2 / j . diameter are 
needed to bring the visual, polarimetric, and ultraviolet photometric estimates 
of the surface pressure of Mars into agreement. The most dramatic evidence 
for the existence of particulate matter in the atmosphere of Mars is the photo- 
graph of the limb of the planet taken by the TV equipment aboard Mariner IV , 
Leighton et al (1965) and Leighton (1966). A cloud estimated to extend to an 
altitude of 100 km was visible in one of the photographs. Nicks (1967) has 
suggested that the cloud might be composed of dry ice crystals . A similar 
suggestion was put forward by Hess (1950) in order to explain the haze layer. 
The F-2 model of Fjeldbo et al (1966a, b) has mesopheric temperatures which 
fall in the range where dry ice crystals would be formed. 

The attachment of electrons and ions to micron and submicron sized 
particles has been observed in the laboratory and the terrestrial atmosphere. 
Formulae for the attachment of charged particles and the resulting charge on 
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the dust particles have been derived by several authors, Fuchs (1947), Gunn 
(1954), Bricard (1962, 1965), and Whipple (1965). 

One of the elements involved in thunderstorm formation is the charging of 
cloud particles see e.g. Shishkin (1965). Such charging should also occur within 
the clouds observed on Mars. The net effect of particulate matter in the Martian 
atmosphere will be a reduction in the electron and ion density distribution. The 
existence of such particle densities at altitudes 20 km above the cloud top is a 
possibility. Attachment to dust should be considered as a potential loss process 
for ions. 

DUSTSTORMS 

It is noticed that yellow clouds move frequently in the vicinity of the thermal 
equator of Mars. Their duration is usually only one to three days and for this 
reason they are termed "duststorms". Newbauer (1966) places the origin of the 
clouds in thermal convection which leads to the formation of dust devils. The 
speed of the wind necessary to initiate the storm has been estimated to be 
145 km/hr and may exceed 300 km/hr, 100 km/hr is frequently observed 
Gifford (1964). Such storms are of interest in the study of atmospheric elec- 
tricity on the planet. For instance Chalmers (1965) has noted that terrestrial 
duststorms give rise to large electric potential gradients and has noted a case 
where lightening was observed. 

CONCLUSION AND RECOMMENDED EXPERIMENTS 

Recent data on the atmosphere and ionosphere of Mars has been utilized 
to outline the electrical properties of the lower Martian atmosphere. Suggested 
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sources of ionization include; solar ultraviolet and x-radiation, cosmic rays, 
meteors, radioactive material within the atmosphere and radioactivity of the 
surface, the solar wind, and protops originating from solar flares. 

Loss processes including charge rearrangement reactions, dissociative 
recombination of ions and electrons and negative ion formation have been de- 
tailed. The principal negative ion may be CO^ if sufficient 0 2 is present. 

The negative ion CO" is also a possibility. 

The altitude distribution of the cosmic ray produced ion density distribution 
is shown to decrease from a maximum of 3 x 10 3 cm -3 at the surface to less 
than 10 2 cm -3 between 55 and 90 km depending on atmospheric model. The 
tail of the solar photoionization rate when combined with this cosmic ray 
production will lead to the formation of an ion density minimum . The magnitude 
of the density and location of the minimum will depend on the atmospheric model 
chosen, the activity of the sun and solar zenith angle as well as such additional 
causes as the presence of particulate matter and meteor ionization. This 
minimum will disappear for Fregian models if solar Lyman alpha plays any 
role in the ionization of the Martian atmosphere. Such an effect could occur if 
nitric oxide is a constituent of the Martian mesosphere. If no negative ions are 
present , the electron density distribution will closely resemble the positive ion 
density curves given in figure 4. The electron density will be less than 5 x 10 3 
cm* 3 at altitudes below 50 km. 

The available data on the Martian ionosphere has lead to several different 
interpretations of the formation of ionosphere, none of which is without objec- 
tions. It is clear that additional experiments will have to be performed before 
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the uncertainty can be resolved. It has been suggested that additional data be 
obtained from an orbiting vehicle. Performance of the occultation experiment 
at smaller solar zenith angles would certainly be useful since under conditions 
the ionization at lower altitudes will be enhanced allowing one to obtain more in- 
formation at lower altitudes. It is such low altitude data which holds the 
solution to the problem of what type ionospheric layer was observed by 
Mariner IV. In order to investigate the lower ionosphere of Mars it is neces- 
sary to launch probes into the atmosphere. Such probes should be instrumented 
to measure concentrations of positive and negative ions as well as their species. 
An experiment should be carried out to monitor the conductivity of the Martian 
atmosphere near the surface of the planet on an extended basis. 

Recent planetary probes to Venus have determined that this planet's 
atmosphere is also largely C0 2 so that the results of the present study can 
be carried over to that atmosphere taking into account, of course the greater 
pressures and higher temperatures encountered on Venus. Vinogradov et al 
(1968) have indicated that the Venusian atmosphere contains water vapor. The 
presence of sufficient quantities of water molecules in the Martian and Venusian 
atmospheres can change the ion species by forming hydrated molecules such 
as C0 4 • (H 2 0)‘ and positive ions of the type H 5 0 2 + and H 3 0 + . 


FIGURE CAPTIONS 


Figure 1 . Electron concentration vs height profile observed with two frequency 
occultation experiment aboard Mariner IV . 

Figure 2. Ion pair production function resulting from solar ultraviolet and 
x-radiation and galactic cosmic radiation E and F-l solar production functions 
are from McElroy (1967). 

Figure 3. Negative ion to electron ratios for different atmospheric models. 

Figure 4. The lower ionosphere of Mars for a) E -region model and b) F-2 
region model. 



REFERENCES 


Alexander, W. M., C. W. McCracken and J. L. Bohn — Science, 149 , 1240-1241, 
(1965) 

Belton, M. J. S. and D. M. Hunten — Astrophys J. 145 , 454-467, (1966) 

Berkner, L. V. and L. C. Marshall — Advances in Geophysics 1967 Academic 
Press, New York 

Bricard, J. — Geofis Pura Appl. 5 1, 237-242, (1962) 

Bricard, J., 82-117 — Problems of Atmospheric and Space Electricity 
Elsevier, Amsterdam (1965) 

Burch, D. S., S, J. Smith and L. M. Brapscomb — Phys. Rev. 112 , 171 (1958) 
Chalmers, J. A., 167-173 — Atmospheric and Space Electricity , Elsevier, 
Amsterdam, (1965) 

Chamberlain, J. W. and M. B. McElroy — Science 152 , 21-25, (1966) 

Chanin, L. M., A. V. Phelps, and M. A. Biondi - Phys. Rev. 128 , 219-230, (1962) 
Donahue, T. M. — Science, 152 , 763-76^, (1966) 

Fehsenfeld, F. C., E. E. Ferguson and A. L. Schmeltekopf — Bui. Am. Phys . 

Soc. Series II 13, 2012, (1968) 

Ferguson, E. E., F. C. Fehsenfeld, P. D. Goldan, A. L. Schneltekopf and H. I. 

Schiff — Planetary Space Sci. 13 , 823-827, (1965) 

Fiocco, G. — J. Geophys. Res., 72 , 3497-3501, (1967) 

Fjeldbo, G., W. C. Fjeldbo and Von R. Eshleman — Science 153 , 1518-1523, 
(1966a) 

Fjeldbo, G., W. C. Fjeldbo and Von R. Eshleman, J. Geophys. Res. 71 , 2307- 
2316, (1966b) 

Fuchs, N. A. — U. S. Army Chemical Laboratories Special Publication CWL 
4:12 (1955) 


REFERENCES (Continued) 


Gray, L. D. — Icarus, v, 390-398, (1966) 

Gifford, F. A. — Monthly Weather Reviews 92, 435-440, (1964) 

Greenspan, J. A. — Science, 150 , 1156-1158, (1965) 

Gunn, R. - J. Meteorol. 11 , 339-347, (1954) 

Hess, S. L. — J. Meteorol. 7 , 1-13, (1950) 

Hintereger, H. E., L. A. Hall and G. Schmidtke — Space Research V , North- 
Holland PublishingCo., Amsterdam (1965) 

Israel, H. — AtmosphSrische ElektrizitSt Akademie Verlag, Leipzig (1961) 
Johnson, F. S. — Science, 150, 1445-1448, (1965) 

Kliore, A., D. L. Cain, G. S. Levy, Von R. Eshleman, G. Fjeldbo, F. D. Drake, 
Science, 149, 1243-1248, (1965) 

Lazarus, A. J., H. S. Bridge, J. M. Davis and C. W. Snyder — 1296-1305 Space 
Research VII, North-Holland Publishing Co., Amsterdam (1967) 

Leighton, R. B., B. C. Murray, R. P. Sharp, J. D. Allen and R. H. Sloan — 
Science 149 , 627-630, (1965) 

Leighton, R. B. — Sci. Am. 214 , 54-68, April (1966) 

Lovell, B. A. — Handbuch der Physik 48 , 427-454, (1957) 

Mariner Stanford Group — Science, 158, 1678-1683, (1967) 

Me Elroy, M. B. - Astrophys J. 150, 11125-1138, (1967) 

McKinley, D. W. R. — Meteor Science and Engineering, McGraw Hill, New York 

J 

(1961) 

Mehr, F. J. and M. A. Biondi — Bui. Am. Phys. Soc. Series II 13, 199, (1968) 
Mentzoni, M. H. and J. Donohoe — Bui. Am. Phys. Soc. Series II 13, 199, (1968) 


REFERENCES (Continued) 


Moruzzi, J. L. and A. V. Phelps — J. Chem. Phys. 45 , 4617-4627, (1966) 
Narcisi, R. S. — Space Research VII I, North-Holland Publishing Co., 
Amsterdam, (1967a) 

Narcisi, R. S., A. D. Bailey and L. Della Lucca — Space Research VII, North- 
Holland Publishing Co., Amsterdam, (1967b) 

Neher, H. V. and H. R. Anderson — J. Geophys. Res. 69 , 807-814, (1964) 
Neubauer, F. M. — J. Geophys. Res. 71, 2419-2426, (1966) 

Nicks, O. W. — A Review of the Mariner IV Results National Aeronautics and 
Space Administration SP-130 (1967) 

Nicolet, M. — Meteors Pergamon Press, London (1955) 

Norton, R. B., E. E. Ferguson, F. C. Fehsenfeld and A. L. Schmeltekopf — 
Planet. Space Sci. 14, 969-979, (1966) 

O’Gallagher, J. J. and J. A. Simpson — Science 149 , 1233-1239, (1965) 

Opik, E. J. — Science, 153 , 255-265, (1966) 

Parker, E. N. - SCI. Am., 210 , 66-80, April (1964) 

Pack, J. L., R. E. Voshall, and A. V. Phelps - Phys. Rev. 127 , 2084-2089, 

(1962) 

Pack, J. L. and A. V. Phelps - J. Chem. Phys. 44 , 1870-1883, (1966a) 

Pack, J. L. and A. V. Phelps — J. Chem. Phys. 42 , 4316-4329, (1966b) 

Sagan, C. and J. Veverka — Science, 158 , 110-112, (1967) 

Schultz, E. D. and Holland Geophysics Corporation of America Technical Report 


62-14-N (1962) 


REFERENCES (Continued) 


Schultz, G. Z. — Phys. Rev. 128 , 178-186 (1962) 

Schlumbohm, H. — Z. Physik, 166 , 192-206, (1962) 

Shishkin, N. S. — 268-279, Problems of Atmospheric and Space Electricity , 
Elsevier, Amsterdam, (1965) 

Smith, E. J., L. Davis, Jr., P. J. Coleman, Jr. — Science, 149 , 1241-1242, 
(1965) 

Van Allen, J. A., L. A. Frank, S. M. Krimigis and H. K. Hills, Science, 149 , 
1228-1233, (1965) 

Vinogradov, A. P., Yu. A. Surkov, K. P. Florenskiy and B. M. Andreychikov 
Doklady A. N. SSSR Astronomiya 179 , No. 1, 37-40, (1968) 

Wameck, P. — J. Chem. Phys. 46, 513-519 (1967a) 

Warneck, P. — Planet. Space Sci, 5 , 1349-1359 (1967b) 

Weller, C. S. and M. A. Biondi, Phys. Rev. Letters 19 , 59-61, (1967) 

Weller, C. S. and M. A. Biondi, Bui. Am. Phys. Sci . Series II, 13, 199, (1968) 
Whipple, E. C. — 123-140, Problems of Atmospheric and Space Electricity 


Elsevier, Amsterdam (1965) 




in 

O 



■wm Banunv 



250 



ION PRODUCTION RATE (cm-3 S ec-i) 





NEGATIVE ION/ELECTRON DENSITY (X) 







CO 


wx~3aruinv 


NUMBER DENSITY* cm' 




wvi-Baniinv 


NUMBER DENSITY- cm 1- 



